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A  new  configuration  of  combined  thermoelectric  device,  two-stage  thermoelectric  refrigerator  driven  by 
two-stage  thermoelectric  generator,  is  proposed  in  this  paper.  The  thermodynamic  model  of  the  com¬ 
bined  device  is  built  by  using  non-equilibrium  thermodynamic  theory.  The  analytical  formulae  for  the 
stable  working  electrical  current,  the  cooling  load  versus  the  working  electrical  current,  and  the  coeffi¬ 
cient  of  performance  (COP)  versus  the  working  electrical  current  of  the  combined  device  are  derived. 
For  the  fixed  total  number  of  thermoelectric  elements  of  the  combined  device,  the  allocations  of  the  ther¬ 
moelectric  element  pairs  among  the  two  thermoelectric  generators  and  the  two  thermoelectric  refriger¬ 
ators  are  optimized  for  maximum  cooling  load  and  COP,  respectively.  The  influences  of  the  heat  source 
temperature  of  the  two-stage  thermoelectric  generator  and  the  heat  source  (cooling  space)  temperature 
of  the  two-stage  thermoelectric  refrigerator  on  the  optimal  performance  of  the  combined  thermoelectric 
device  are  analyzed  by  detailed  numerical  examples. 
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1.  Introduction 

Semiconductor  thermoelectric  power  generation,  based  on  the 
Seebeck  effect,  and  semiconductor  thermoelectric  cooling,  based 
on  the  Peltier  effect,  have  very  interesting  capabilities  compared 
to  conventional  power  generation  and  cooling  systems  [1-4].  The 
absence  of  moving  components  results  in  an  increase  of  reliability, 
a  reduction  of  maintenance,  and  an  increase  of  system  life;  the 
modularity  allows  for  application  in  a  wide-scale  range  without 
significant  losses  in  performance;  the  absence  of  a  working  fluid 
avoids  environmental  dangerous  leakage;  and  the  noise  reduction 
appears  also  to  be  an  important  feature.  Thermoelectric  generator 
and  refrigerator  have  been  used  in  military,  aerospace,  instrument, 
and  industrial  or  commercial  products,  as  a  power  generation  and 
cooling  devices  for  specific  purposes.  Many  researchers  are  con¬ 
cerned  about  the  physical  properties  of  thermoelectric  material 
and  the  manufacturing  technique  of  thermoelectric  modules  [5- 
8].  In  addition  to  the  improvement  of  the  thermoelectric  material 
and  module,  the  system  analysis  and  optimization  of  thermoelec¬ 
tric  generator  and  refrigerator  are  equally  important  in  designing 
high-performance  thermoelectric  generators  and  refrigerators. 

In  general,  conventional  non-equilibrium  thermodynamics 
[1,9,10]  is  used  to  analyze  the  performance  of  single-stage  one- 
or  multiple-element  thermoelectric  generators  [11-20]  and  refrig¬ 
erators  [21-32].  Due  to  the  performance  limits  of  thermoelectric 
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material,  thermoelectric  generators  and  refrigerators  of  two  stages 
or  more  should  be  applied  to  improve  the  level  of  thermodynamic 
performance.  The  performance  analysis  and  optimization  of  two- 
stage  thermoelectric  generators  [33,34]  and  refrigerators  [35-41] 
were  also  performed.  All  of  those  were  performed  by  using  conven¬ 
tional  non-equilibrium  thermodynamics  without  considering  the 
external  losses.  Performances  of  two-stage  thermoelectric  genera¬ 
tors  [42]  and  refrigerators  [43]  with  external  heat  transfer  are  ana¬ 
lyzed  using  the  combination  of  finite  time  thermodynamics  and 
non-equilibrium  thermodynamics  by  Chen  et  al.  [42,43]. 

All  objects  of  these  researches  are  dependent  thermoelectric  de¬ 
vices,  that  is,  they  have  a  direct-current  power  source  providing  di¬ 
rect  current  to  refrigerate  or  heat  up.  However,  for  some  special 
systems,  such  as  submarines,  cars,  and  special  electric  circuit,  the 
heat  rejected  from  the  thermal  machine  may  drive  a  thermoelec¬ 
tric  refrigerator  through  the  use  of  a  thermoelectric  generator,  so 
that  the  thermoelectric  cooler  does  not  need  an  independent 
power  source.  Such  a  new  refrigeration  system  is  directly  com¬ 
posed  of  a  thermoelectric  generator  and  a  thermoelectric  cooler. 
It  is  different  from  the  traditional  thermoelectric  systems  which 
merely  consist  of  thermoelectric  generators  or  coolers,  and  de¬ 
serves  to  be  investigated  both  from  the  point  of  view  of  theoretical 
design,  as  well  as  practical  application.  Chen  et  al.  [44]  and  Khattab 
and  El  Shenawy  [45]  built  a  model  of  this  kind  of  combined  system, 
single-stage  thermoelectric  refrigerator  driven  by  single-stage 
thermoelectric  generator,  and  analyzed  the  performance  of  the  de¬ 
vice.  Based  on  the  performance  analysis  and  optimization  of  ther¬ 
moelectric  generator  and  heat  pump  by  using  non-equilibrium 
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thermodynamics,  Meng  et  al.  [46]  built  a  model  of  single-stage 
thermoelectric  heat  pump  driven  by  single-stage  thermoelectric 
generator,  and  analyzed  and  optimized  its  performance. 

Based  on  the  performance  analysis  and  optimization  of  two- 
stage  thermoelectric  generators  [33,34]  and  two-stage  thermoelec¬ 
tric  refrigerators  [35-41]  by  using  non-equilibrium  thermodynam¬ 
ics,  this  paper  provides  a  new  configuration  of  combined 
thermoelectric  device,  two-stage  thermoelectric  refrigerator  dri¬ 
ven  by  two-stage  thermoelectric  generator.  The  thermodynamic 
model  of  the  combined  device  is  built  by  using  non-equilibrium 
thermodynamic  theory.  Three  analytical  formulae  for  the  stable 
working  electrical  current,  the  cooling  load  versus  the  working 
electrical  current,  and  the  coefficient  of  performance  (COP)  versus 
the  working  electrical  current  of  the  combined  device  are  derived. 
For  a  fixed  total  number  of  thermoelectric  elements  of  the  com¬ 
bined  device,  the  thermoelectric  element  allocations  among  the 
two  thermoelectric  generators  and  the  two  thermoelectric  refriger¬ 
ators  are  optimized  for  maximizing  cooling  load  and  COP, 
respectively.  There  has  been  no  investigation  concerning  the  per¬ 
formance  analysis  and  optimization  for  such  combined  thermo¬ 
electric  devices  published  in  the  open  literature.  The  influences 
of  the  heat  source  temperature  of  the  two-stage  thermoelectric 
generator  and  the  heat  source  (cooling  space)  temperature  of  the 
two-stage  thermoelectric  refrigerator  on  the  optimal  performance 
of  the  combined  device  are  analyzed  by  detailed  numerical 
examples. 


couple.  The  conduction  heat  losses  are  I<{Tm  -  Tm )  for  the  top  stage 
and  I<(Tm  -  TL i)  for  the  bottom  stage,  respectively,  where  I<  is  the 
thermal  conductance  of  the  semiconductor  couple,  Tm  is  the  hot 
junction  (heat  source)  temperature,  TL i  is  the  cold  junction  (heat 
sink)  temperature,  and  Tm  is  the  temperature  of  the  junction  be¬ 
tween  the  two  stages.  The  rates  of  heat  flow  of  the  thermoelectric 
generator  are  Qhi,  Qm,  and  On. 

The  refrigerator  consists  of  a  top  stage  with  nA  pairs  of  thermo¬ 
electric  elements  and  a  bottom  stage  with  n2  pairs  of  thermoelec¬ 
tric  elements.  The  total  number  of  thermoelectric  element  pairs  of 
the  refrigerator  is  n,  i.e.  n  =  +  n2.  Each  element  is  composed  of  a 

P-type  and  a  N-type  semiconductor  legs.  The  structure  of  the 
refrigerator  is  similar  to  the  generator.  The  conduction  heat  losses 
are  I<{TH 2  -  Tn )  for  the  top  stage  and  I<{Tn  -  TL2 )  for  the  bottom 
stage,  respectively,  where  TH2  is  the  hot  junction  (heat  sink)  tem¬ 
perature,  TL2  is  the  cold  junction  (heat  source)  temperature,  and 
Tn  is  the  temperature  of  the  junction  between  the  two  stages. 
The  rates  of  heat  flow  of  the  thermoelectric  refrigerator  are  Qh2, 
Qn,  and  Ql2. 

The  total  number  of  thermoelectric  element  pairs  of  the  com¬ 
bined  thermoelectric  device,  M,  is  finite  and  M  =  m  +  n  holds.  The 
cooling  load  of  the  combined  thermoelectric  device  is  Ql2-  The 
power  input  required  by  the  refrigerator  is  the  power  output  of 
the  generator. 

According  to  the  theory  of  non-equilibrium  thermodynamics, 
for  the  two-stage  thermoelectric  generator,  one  has 


2.  Model  of  a  combined  thermoelectric  device 

A  schematic  diagram  of  a  combined  thermoelectric  device  is 
shown  in  Fig.  1.  The  device  consists  of  a  two-stage  thermoelectric 
generator  and  a  two-stage  thermoelectric  refrigerator  in  series.  The 
direct-current  power  source  of  the  refrigerator  is  the  current  of  the 
generator. 

The  generator  consists  of  a  top  stage  with  mi  pairs  of  thermo¬ 
electric  elements  and  a  bottom  stage  with  m2  pairs  of  thermoelec¬ 
tric  elements.  The  total  number  of  thermoelectric  element  pairs  of 
the  generator  is  m,  i.e.  m  =  m\  +  m2.  Each  element  is  composed  of  a 
P-type  and  a  N-type  semiconductor  legs.  The  thermoelectric  gener¬ 
ating  element  is  assumed  to  be  insulated,  both  electrically  and 
thermally,  from  its  surroundings,  except  at  the  junction-reservoir 
contacts  and  the  junction  between  the  two  stages.  The  internal 
irreversibility  is  caused  by  Joule  loss  and  heat  conduction  loss 
through  the  semiconductor  between  the  hot  and  cold  junctions. 
The  Joule  loss  generates  an  internal  heat  I2R,  where  R  is  the  total 
internal  electrical  resistance  of  the  semiconductor  couple  and  I  is 
the  working  electrical  current  generating  from  the  semiconductor 


Qhi  —  m\ 

Qm=m  i 
Qm=™2 

Ql  i  =  m2 


aITm+K(Tm-Tm)-h2R 


<xITm  +  I((Tm  —Tm)+—rR 
ocITm+K(Tm-TL1)-^I2R 
cdTn+K(Tm-TL,)+h2R 


For  the  two-stage  thermoelectric  refrigerator,  one  has 


a) 

(2) 

(3) 

(4) 


0.H2  ~  nl 
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(5) 
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(6) 
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Fig.  1.  Schematic  diagram  of  the  combined  thermoelectric  device. 
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where  a  =  aP  -  ccN,  aP,  and  ocN  are  the  Seebeck  coefficients  of  the  P- 
and  N-type  semiconductor  legs  for  each  thermoelectric  power  gen¬ 
eration  or  refrigerating  element. 

Combining  Eq.  (2)  with  Eq.  (3)  gives  Tm.  Substituting  it  into  Eqs. 
(1)  and  (4)  yields  Qhi  and  Q± i-  Combining  Eq.  (6)  with  Eq.  (7)  gives 
Tn.  Substituting  it  into  Eqs.  (5)  and  (8)  yields  Qh2  and  Ql2.  According 
to  energy  balance  equation  Qm  +  Ql2  =  On  +  Qh2,  one  can  obtain  Eq. 
(A8)  that  should  be  satisfied  by  a  stable  electrical  current,  see 
Appendix  A. 

For  the  given  parameters  a,  K,  R,  Tm,  Tu,  TH2,  TL2,  m^,  m2 ,  Hi,  n2, 
one  can  obtain  the  system  stable  current  Is  (see  Appendix  A),  and 
then  obtain  the  cooling  load  and  COP  of  the  combined  thermoelec¬ 
tric  device  as  follows: 

0.L2  =  n2{ocIsTL2  -  JC[(0.5ni/s2R  +  0.5 n,/2R  +  n^KTm 

+  n2KTH2)/(n1a/s  -  n2a/s  +  n^K  +  n2K)  -  TL2]  -  0.5 12SR}  (9) 


sQl2.  In  general,  the  optimum  allocations  of  thermoelectric  element 
pairs  and  optimum  working  currents  at  maximum  cooling  load  are 
different  from  these  at  maximum  COP.  The  maximum  cooling  load 
of  the  combined  thermoelectric  device,  Qi2imax,  is  directly  propor¬ 
tional  to  the  total  number  of  thermoelectric  element  pairs  (see 
Fig.  2),  and  the  maximum  COP,  emax,  is  independent  of  the  total 
number  of  thermoelectric  element  pairs.  The  optimum  working 
currents  corresponding  to  the  maximum  cooling  load,  Iopt,aL2,  and 
the  optimum  working  currents  corresponding  to  the  maximum 
COP,  /opt,e,  are  also  independent  of  the  total  number  of  thermoelec¬ 
tric  element  pairs,  and  Iopt,QL2  >  /opt,e  holds.  The  cooling  load  at 
maximum  COP,  Qz.2e,  increases  with  the  increase  of  the  total  num¬ 
ber  of  thermoelectric  element  pairs,  while  the  COP  at  maximum 
cooling  load,  £Ql2,  is  independent  of  the  total  number  of  thermo¬ 
electric  element  pairs. 

Although  the  optimum  thermoelectric  element  allocations 
among  the  two  thermoelectric  generators  and  the  two  thermoelec- 


8  —  Q.L2/Q.H1 


n2{«/srL2  -  K[(0.5ni  J2R  +  O^n^R  +  ihKTm  +  n2KTH2)/(nia/s  -  n2a/s  +  n^K  +  n2/Q  -  Tq]  -  0.5/2R} 

mi{a/srm  +  fC[rH1  -  (0.5m,/2R  +  0.5m2/2R  +  m,/<TH1  +  m2IO'n)/(m2aIs  -  m,a/s  +  m,fC  +  m2K)}  -  0.5I2SR} 


(10) 


3.  Numerical  examples 

Numerical  calculations  are  performed  in  order  to  analyze  and 
optimize  the  performance  of  the  two-stage  thermoelectric  refriger¬ 
ator  driven  by  two-stage  thermoelectric  generator.  The  total  num¬ 
ber  of  thermoelectric  element  pairs  is  finite,  i.e.  M  =  m1+m2  + 
xi\  +  n2.  The  optimum  thermoelectric  element  allocations  and  the 
optimum  working  electrical  current,  Iop t,aL2  or  W£  at  maximum 
cooling  load,  QL2,max,  or  maximum  COP,  emax,  are  searched,  respec¬ 
tively.  In  the  calculations,  a  =  2.1  x  1(T4V/K,  I<=  1.6  x  10“2  W/K, 
and  R  =  1.2  x  10-3  Q  are  used  [47]. 

3.1.  Effects  of  total  number  and  allocation  of  thermoelectric  element 
pairs 

Table  1  lists  the  optimum  parameters  and  device  performance 
with  different  total  number  of  thermoelectric  element  pairs.  In 
the  calculations,  the  temperatures  of  the  heat  reservoirs  are  set 
to  be  Tm  =  450  K,  TL1  =  300  K,  TH2  =  300  K,  and  TL2  =  280  K.  One 
can  see  from  the  table  that  there  exist  optimum  thermoelectric  ele¬ 
ment  allocations  among  the  two  thermoelectric  generators  and  the 
two  thermoelectric  refrigerators  for  the  fixed  total  number  of  ther¬ 
moelectric  element  pairs  and  there  exist  optimum  working  cur¬ 
rents,  Iopt,dL2  or  J0pt,e,  corresponding  to  maximum  cooling  load, 
0L2,max*  or  maximum  COP,  emax,  respectively.  The  cooling  load  at 
maximum  COP  is  Qi2.e,  and  the  COP  at  maximum  cooling  load  is 


trie  refrigerators  are  different  for  different  total  number  of  thermo¬ 
electric  element  pairs  and  different  optimizing  objectives,  there  are 
some  general  rules.  Three  ratios  of  number  of  thermoelectric  ele¬ 
ment  pairs  are  defined:  total  element  ratio  of  generator  to  refriger¬ 
ator,  x  =  m/M,  generator  element  ratio,  Xi  =  mi/m,  and  refrigerator 
element  ratio,  x2  =  ni/n.  Table  2  lists  the  results  of  optimum  ratios 
of  number  of  thermoelectric  element  pairs.  One  can  see  that  three 
optimum  ratios  of  number  of  thermoelectric  element  pairs  for 
maximum  cooling  load,  xopt,Qi2,  *ioPt,QL2,  and  x2opt,QL2,  and  three  opti¬ 
mum  ratios  of  number  of  thermoelectric  element  pairs  for  maxi¬ 
mum  COP,  xopt,E,  *iopt.£»  and  x2optl£,  are  independent  of  the  total 
number  of  thermoelectric  element  pairs.  If  mi,  m2,  nlt  and  n2  in 
Eqs.  (A9)-(A12)  are  replaced  by  M,  x,  xlf  and  x2,  one  can  see  that 
the  equation  for  the  stable  working  electrical  current,  Eq.  (A8)  is 
independent  of  the  total  number  of  thermoelectric  element  pairs, 
M.  In  theory,  given  x,  Xi,  and  x2,  the  total  number  of  thermoelectric 
element  pairs,  M  does  not  affect  the  working  electrical  current.  In 
practice,  the  numbers  of  the  thermoelectric  element  pairs  of  the 
two  generators  and  the  two  refrigerators,  mlt  m2,  ni,  and  n2  must 
be  integers,  which  makes  the  optimum  ratios,  xopt,  Xiopt,  and  x2opt 
vary  a  little  with  the  total  number  of  thermoelectric  element  pairs, 
M.  In  this  example,  one  can  see  that  xopt,QL2  «  0.64  and  xopt,e  «  0.51. 
That  is,  xopt,Qi2  >  x0pt,£,  the  optimum  total  element  ratio  of  generator 
to  refrigerator  at  maximum  cooling  load  is  larger  than  that  at  max¬ 
imum  COP.  This  means  that  the  optimum  allocation  range  is 
*oPt,£  ^x^xopt,QL2,  more  thermoelectric  elements  in  the  thermo- 


Table  1 

Optimum  parameters  and  performance  of  the  combined  thermoelectric  device. 
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Fig.  2.  Maximum  cooling  load  versus  total  number  of  thermoelectric  element  pairs. 


electric  generators  can  obtain  larger  cooling  load,  and  more  ther¬ 
moelectric  elements  in  the  thermoelectric  refrigerators  can  obtain 
larger  COP.  For  the  optimum  total  element  ratio  of  generator  to 
refrigerator,  either  for  maximum  cooling  load  or  for  maximum 
COP,  the  optimum  generator  element  ratios  are  always  Xiopt,Qi2  = 
XioPt,e  =  0.5,  i.e.  the  equal  partition  of  number  of  thermoelectric 
element  between  the  two  thermoelectric  generators  will  lead  to 
both  maximum  cooling  load  and  maximum  COP.  While  the 
optimum  refrigerator  element  ratios,  x2opt,Qi2  and  x2opt,£  are  alwaYs 
larger  than  0.5,  they  are  also  almost  the  same. 

3.2.  Effects  of  heat  source  temperature  of  the  two-stage  thermoelectric 
generator 

The  effects  of  the  heat  source  temperature  of  the  two-stage 
thermoelectric  generator  on  the  optimum  parameters  and  the  opti¬ 
mum  performance  of  the  combined  thermoelectric  device  are  ana¬ 
lyzed.  In  the  calculations,  Tn  =  300  K,  TH2  =  300  K,  TL2  =  280  K,  and 
M  =  200  are  set.  Figs.  3  and  4  show  the  characteristic  that  three 
optimum  ratios  of  number  of  thermoelectric  element  pairs  for 
maximum  cooling  load  (xopt,Qi2 ,  Xiopt,Qi2 ,  and  x2opt,Qi2 )  versus  the 
heat  source  temperature  of  the  two-stage  thermoelectric  generator 
and  three  optimum  ratios  of  number  of  thermoelectric  element 
pairs  for  maximum  COP  (xopt,e,  Xiopt,e,  and  x2opt,e)  versus  the  heat 
source  temperature  of  the  two-stage  thermoelectric  generator, 
respectively.  The  corresponding  maximum  cooling  load  and  maxi¬ 
mum  COP  versus  the  heat  source  temperature  of  the  two-stage 
thermoelectric  generator  are  shown  by  solid  lines  in  Figs.  5  and 
6,  respectively.  For  comparisons,  the  cooling  load  of  the  combined 


Fig.  3.  Optimum  parameters  at  maximum  cooling  load  versus  the  heat  source 
temperature  of  the  two-stage  generator. 


thi  (K) 


Fig.  4.  Optimum  parameters  at  maximum  COP  versus  the  heat  source  temperature 
of  the  two-stage  generator. 


thermoelectric  device  with  mA  =  53,  m2  =  53,  n\  =  56,  and  n2  =  38, 
which  are  the  optimum  numbers  of  the  thermoelectric  element 
pairs  of  the  combined  thermoelectric  device  for  maximum  cooling 
load  objective  with  Tm  =  600  K,  is  shown  by  dotted  line  in  Fig.  5, 


Table  2 

Optimum  ratios  of  number  of  thermoelectric  element  pairs. 
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Fig.  5.  Maximum  cooling  load  versus  the  heat  source  temperature  of  the  two-stage 
generator. 


tHi  (X> 


Fig.  6.  Maximum  COP  versus  the  heat  source  temperature  of  the  two-stage 
generator. 

and  the  COP  of  the  combined  thermoelectric  device  with  mi  =  33, 
m2  =  33,  Hi  =  76,  and  n2  =  58,  which  are  the  optimum  numbers  of 
the  thermoelectric  element  pairs  of  the  combined  thermoelectric 
device  for  maximum  COP  objective  with  Tm  =  600  K,  is  shown  by 
dotted  line  in  Fig.  6.  The  corresponding  optimum  working  electri¬ 
cal  currents,  JQi2,  for  maximum  cooling  load  and  Je,  for  maximum 
COP  versus  the  heat  source  temperature  of  the  two-stage  thermo¬ 
electric  generator  are  shown  in  Fig.  7. 

One  can  see  that  the  optimum  total  element  ratio  of  generator 
to  refrigerator  xopt,Qi2  and  xopt,£  decrease  with  the  increase  in  the 
heat  source  temperature  of  the  two-stage  thermoelectric  genera¬ 
tor.  This  means  that  more  thermoelectric  element  pairs  must  be 
allocated  to  the  refrigerator  in  order  to  obtain  maximum  cooling 
load  or  maximum  COP  when  the  heat  source  temperature  of  the 
two-stage  thermoelectric  generator  increases.  The  selection  range 
of  the  total  element  ratio  of  generator  to  refrigerator  should  be 
x8  ^x^xQi2.  More  thermoelectric  element  pairs  in  the  thermo¬ 
electric  generator  can  obtain  larger  cooling  load,  and  more  thermo¬ 
electric  element  pairs  in  the  thermoelectric  refrigerator  can  obtain 
larger  COP. 

The  optimum  generator  element  ratio  at  maximum  cooling 
load,  Xiopt,Qi2  and  the  optimum  generator  element  ratio  at  maxi¬ 
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Fig.  7.  Optimum  working  currents  versus  the  heat  source  temperature  of  the  two- 
stage  generator. 


mum  COP,  XioPt,e,  are  always  0.5  for  different  heat  source  tempera¬ 
ture  of  the  two-stage  thermoelectric  generator. 

The  optimum  refrigerator  element  ratios,  x2opt,Qi2  and  x2opt,£, 
change  little  with  the  increase  of  the  heat  source  temperature  of 
the  two-stage  thermoelectric  generator. 

Both  the  maximum  cooling  load  and  the  maximum  COP  in¬ 
crease  monotonically  with  the  increase  of  the  heat  source  temper¬ 
ature  of  the  two-stage  thermoelectric  generator  (See  solid  lines  in 
Figs.  5  and  6).  If  the  allocations  of  the  thermoelectric  element  pairs 
among  the  two  thermoelectric  generators  and  the  two  thermoelec¬ 
tric  refrigerators  do  not  meet  the  needs  of  optimum  allocations  for 
maximum  cooling  load  or  maximum  COP,  both  cooling  load  and 
COP  are  less  than  the  maximum  ones. 

Both  the  optimum  working  currents  IqL2  and  IE  increase  with  the 
increase  of  the  heat  source  temperature  of  the  two-stage  thermo¬ 
electric  generator,  and  IQl2  ^  Is  holds.  The  difference  between  7Qi2 
and  IE  increase  with  the  increase  of  Tm  and  if  Tm  -►  Tn  =  300  K, 
IQl2  -  Is  =  Td  — >  0.  The  optimum  working  currents  for  both  heating 
load  and  COP  are  Is  <  Jopt  ^  IqL2. 

3.3.  Effects  of  heat  source  temperature  of  the  two-stage  thermoelectric 
refrigerator 

The  effects  of  the  heat  source  (cooling  space)  temperature  of  the 
two-stage  thermoelectric  refrigerator  on  the  optimum  parameters 
and  the  optimum  performance  of  the  combined  thermoelectric  de¬ 
vice  are  analyzed.  In  the  calculations,  Tm  =  800  K,  TL i  =  300  K, 
Th 2  =  300  K,  and  M  =  200  are  set.  Figs.  8  and  9  show  the  character¬ 
istic  that  three  optimum  ratios  of  number  of  thermoelectric 
element  pairs  for  maximum  cooling  load  (xopt,Qi2,  *ioPt,QL2, 
and  x2opt,Qi2 )  versus  the  heat  source  temperature  of  the  two-stage 
thermoelectric  refrigerator  and  three  optimum  ratios  of  number 
of  thermoelectric  element  pairs  for  maximum  COP  (xopti£,  Xioptfe, 
and  x2opt.£)  versus  the  heat  source  temperature  of  the  two-stage 
thermoelectric  refrigerator,  respectively.  The  corresponding  maxi¬ 
mum  cooling  load  and  maximum  COP  versus  the  heat  source  tem¬ 
perature  of  the  two-stage  thermoelectric  refrigerator  are  shown  by 
solid  lines  in  Figs.  10  and  11,  respectively.  For  comparisons,  the 
cooling  load  of  the  combined  thermoelectric  device  with  mi  =  45, 
m2  =  45,  ni  =  68,  and  n2  =  42,  which  are  the  optimum  numbers  of 
the  thermoelectric  element  pairs  of  the  combined  thermoelectric 
device  for  maximum  cooling  load  objective  with  TL2  =  280  K,  is 
shown  by  dotted  line  in  Fig.  10,  and  the  COP  of  the  combined  ther- 
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Fig.  8.  Optimum  parameters  at  maximum  cooling  load  versus  the  heat  source 
temperature  of  the  two-stage  refrigerator. 


Fig.  11.  Maximum  COP  versus  the  heat  source  temperature  of  the  two-stage 
refrigerator. 


Fig.  9.  Optimum  parameters  at  maximum  COP  versus  the  heat  source  temperature 
of  the  two-stage  refrigerator. 
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Fig.  10.  Maximum  cooling  load  versus  the  heat  source  temperature  of  the  two- 
stage  refrigerator. 


moelectric  device  with  mi  =  20,  m2  =  26,  nA  =  87,  and  n2  =  67,  which 
are  the  optimum  numbers  of  the  thermoelectric  element  pairs  of 
the  combined  thermoelectric  device  for  maximum  COP  objective 
with  TL2  =  280  K,  is  shown  by  dotted  line  in  Fig.  11.  The  correspond¬ 
ing  optimum  working  currents,  JQt2  and  IE  versus  the  heat  source 
temperature  of  the  two-stage  thermoelectric  refrigerator  are 
shown  in  Fig.  12. 

One  can  see  that  the  optimum  total  element  ratio  of  generator 
to  refrigerator  xopt,Qi2  and  xopt,e  decrease  with  the  increase  of  the 
heat  source  temperature  of  the  two-stage  thermoelectric  refriger¬ 
ator.  This  means  more  thermoelectric  element  pairs  must  be  allo¬ 
cated  to  the  refrigerator  in  order  to  obtain  maximum  cooling  load 
or  maximum  COP  when  the  heat  source  temperature  of  the  two- 
stage  thermoelectric  refrigerator  increases. 

Both  the  optimum  generator  element  ratio  Xiopt,QL2  at  maximum 
cooling  load  and  the  optimum  generator  element  ratio  Xiopt£  at 
maximum  COP  are  always  0.5  for  different  heat  source  tempera¬ 
ture  of  the  two-stage  thermoelectric  refrigerator. 

The  optimum  refrigerator  element  ratios  x2opt,Qi2  and  x2opt  e  de¬ 
crease  with  the  increase  of  the  heat  source  temperature  of  the 
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Fig.  12.  Optimum  working  currents  versus  the  heat  source  temperature  of  the  two- 
stage  refrigerator. 
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two-stage  thermoelectric  refrigerator.  This  means  that  more  ther¬ 
moelectric  element  pairs  must  be  allocated  to  the  bottom  stage  of 
the  refrigerator  in  order  to  obtain  maximum  cooling  load  or  max¬ 
imum  COP  when  the  heat  source  temperature  of  the  two-stage 
thermoelectric  refrigerator  increases. 

Both  the  maximum  cooling  load  and  the  maximum  COP  in¬ 
crease  rapidly  with  the  increase  of  the  heat  source  temperature 
of  the  two-stage  thermoelectric  refrigerator.  If  the  allocations  of 
the  thermoelectric  element  pairs  among  the  two  thermoelectric 
generators  and  the  two  thermoelectric  refrigerators  do  not  meet 
the  needs  of  optimum  allocations  for  maximum  cooling  load  or 
maximum  COP,  both  the  cooling  load  and  the  COP  are  less  than 
the  maximum  ones. 

Both  the  optimum  working  currents  JQl2  and  IE  decrease  with  the 
increase  of  the  heat  source  temperature  of  the  two-stage  thermo¬ 
electric  refrigerator,  and  IQl2  ^  Is  holds.  The  difference  between  JQi2 
and  IE  increases  with  the  increase  of  TL2  and  if  Tm  200  K  (that  is 
the  minimum  cooling  temperature  in  this  condition),  7Ql2  -  Is  = 
Td  — >  0.  The  optimum  working  currents  for  both  heating  load  and 
COP  are  I8  <  Jopt  <  IqL2. 

4.  Conclusion 

A  new  configuration  of  combined  thermoelectric  device,  two- 
stage  thermoelectric  refrigerator  driven  by  two-stage  thermoelec¬ 
tric  generator,  is  proposed  in  this  paper.  The  thermodynamic  mod¬ 
el  of  the  combined  device  is  built  by  using  non-equilibrium 
thermodynamic  theory.  Three  analytical  formulae  for  the  stable 
working  electrical  current,  the  cooling  load  versus  the  working 
electrical  current,  and  the  COP  versus  the  working  electrical  cur¬ 
rent  of  the  combined  device  are  derived.  The  performance  optimi¬ 
zation  of  the  combined  thermoelectric  device  is  performed  by 
searching  the  allocations  of  the  thermoelectric  element  pairs 
among  the  two  thermoelectric  generators  and  the  two  thermoelec¬ 
tric  refrigerators.  There  has  been  no  investigation  concerning  the 
performance  analysis  and  optimization  for  such  combined  thermo¬ 
electric  device  published  in  the  open  literature.  The  influences  of 
the  heat  source  temperature  of  the  two-stage  thermoelectric  gen¬ 
erator  and  the  heat  source  (cooling  space)  temperature  of  the  two- 
stage  thermoelectric  refrigerator  on  the  optimal  performance  of 
the  combined  thermoelectric  device  are  analyzed.  All  the  parame¬ 
ters  should  be  considered  in  the  design  and  application  of  practical 
combined  thermoelectric  devices  in  order  to  obtain  the  maximum 
economy  benefit. 

The  results  show  that  there  exist  optimum  thermoelectric  ele¬ 
ment  allocations  among  the  two  thermoelectric  generators  and 
the  two  thermoelectric  refrigerators  for  the  fixed  total  number  of 
thermoelectric  element  pairs  and  there  exist  the  optimum  working 
currents  corresponding  to  maximum  cooling  load  or  maximum 
COP,  respectively.  In  general,  the  optimum  allocations  of  thermo¬ 
electric  element  pairs  and  optimum  working  currents  at  maximum 
cooling  load  are  different  from  these  at  maximum  COP.  The  opti¬ 
mum  working  currents  corresponding  to  the  maximum  cooling 
load  is  larger  than  that  corresponding  to  the  maximum  COP.  The 
cooling  load  at  maximum  COP  increases  with  the  increase  of  the 
total  number  of  thermoelectric  element  pairs,  while  the  COP  at 
maximum  cooling  load  is  independent  of  the  total  number  of  ther¬ 
moelectric  element  pairs.  Three  optimum  ratios  of  number  of  ther¬ 
moelectric  element  pairs  for  maximum  cooling  load  and  three 
optimum  ratios  of  number  of  thermoelectric  element  pairs  for 
maximum  COP  are  also  independent  of  the  total  number  of  ther¬ 
moelectric  element  pairs. 

The  results  obtained  herein  may  provide  guidelines  for 
the  design  and  application  of  practical  combined  thermoelectric 
devices. 
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Appendix  A.  A  stable  electrical  current 

Combining  Eq.  (2)  with  Eq.  (3)  gives 

(mi  -f  m2)I2R/ 2  +  m,KTm  +  m 2KT H2 
m  <x/(m2-m,)+K(m,  +  m2)  1  j 

Combining  Eq.  (6)  with  Eq.  (7)  gives 

j  _  (Hi  +  n2)l2R/2  +  ri\KTm  +  h2/<Th2 

n~  a/(ni  -n2)  +  /C(n1+n2)  1  J 

Substituting  Eq.  (Al)  into  Eqs.  (1)  and  (4)  yields 
Q-hi  —  +  I([Tm  —  (0.5tn^I2R  +  0.5m2I2R  +  nri\KTn\ 

+  m2JCTLi)/(m2aJ  -  mi  a/  +  m^K  +  m2I<)]  -  0.5/2R}  (A3) 

QL1  =  m2{a/TL1  +  K[(0.5mi  J2R  +  0.5  m2I2R  +  mAKTm 
+  m2/CTLi)/(m2a/  -  rriia/  +  m^K  +  m2I<)  -  Tn\ 

+  0.5J2R}  (A4) 

Substituting  Eq.  (A2)  into  Eqs.  (5)  and  (8)  yields 
Q-H2  —  Hi  {ocITh2  —  K[Th2  —  (0.5ni/2R  +  0.5ni/2R  +  n^ICT^ 

+  n2/<TH2)/(nia/  -  n2aJ  +  nxK  +  n2/C)]  +  0.5/2R}  (A5) 

qL2  =  n2{a/TL2  -  JC[(0.5ni/2R  +  0.5ni/2R  +  n^KTm 

+  n2/<TH2)/(nia/  -  n2a/  +  n^K  +  n2I< )  -  TL2]  -  0.5/2R}  (A6) 

The  overall  system  is  a  closed  loop  circuit,  and  heat  flow  of  the 
system  is  in  balance,  one  has 

Qni  +  0-L2  =  0-Ll  +  0-H2  (A7) 

Substituting  Eqs.  (A3)-(A6)  into  Eq.  (A7)  and  re-arranging  the 
results  yields  the  equation  that  should  be  satisfied  by  a  stable  elec¬ 
trical  current  I 

A3I3  +A2I2  +A!/  +  A0  =  0  (A8) 

where 

A3  =  a2R(mi  -  m2)(n2  -  ni)(mi  +  m2  +  n  1  +  n2)  (A9) 

A2  =  a[a2(n,  -  n2)(m,  -  m2)(mJH ,  +  n2TL2  -  nJH2 

-  m2Tu  )  +  0.5RK(4n2n,mi  -njm2  -3m^ni  +  n2m\ 

-  m\n\  +3 m\n2  +  n2m^  -  4n2n1m2  +3n^m! 

+  4m2mAn2 -4m1m2ni -3n2m2)]  (A10) 

Ai  =  2a2K(n^ti2m2TH2  -  3m,m2n^Tm  +  m,m2n2Tm 
-3m2mjn2Tu  +m2m,n,TL1  +m2n1TH1  +m2n2Tu 
+  mjnJu  -  3nin2miTH2  +  n22mjl2  +  n2m^TH2 
+  m2n2TH  1  +  n2nimi TL2  - 3n2n^m2TL2  +  njm2TL2 
+  n\m2TH2)  -  2 RI<2(n,  +  n2)(mt  +  m2)(mt  +m2+  n, 

+  n2)  (All) 
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A0  =  4(xK2(-mxTHXm2n2  +  nxTH2n2mx  +  nxTH2n2m2 

-  mxTHXm2nx  +  m2TLXmxnx  +  m2TLXmxn2 

-  n2Tl2nxmx  -  n2TL2nxm2)  (A12) 

The  device  practical  electrical  current  (the  solution  of  Eq.  (A8))  is 
dependent  on  the  values  of  mx,  m2 ,  nx  and  n2. 

When  mx  =  m2  or  nx  =  n2,  then  one  has  that  A3=A2  =  0  and  the 
practical  electrical  current  Is  is 

4=/°  =  -^  (A13) 


When  mx  ^  m2  and  nx  ^  n2 ,  the  solutions  of  Eq.  (A8)  are 


Ix  =  Bx  +  B2  -\-^-B3i 

(A14) 

h  =  Bx  +  B2  — 

(A15) 

h  =  B2  +  B3 

(A16) 

where 

C1/3  3AxA3-A22 

1  12  a3  +  3a3c1/3 

(A17) 

D  1*2 

B2  =  ~3  A3 

(A18) 

D  C1/3  ,  6A,A3  -  2A3 

B3~6A3+  3A3c^ 

(A19) 

and 

C  =  36 A3A2A3  -  108A0A3  -  8a3  +  nV3A3(4A2A3  -  A\A\ 

-  18AoA,A2A3  +  27A20A23  +  4A0A3)1/2 

(A20) 

Analysis  shows  that  the  practical  stable  electrical  current  solu¬ 
tion  Is  must  satisfy  the  following  rules: 


(1)  When  mi 
h  =  I,  =  B , 

<  m2  and  nx  >  n2, 
+  B2+^B3i. 

the 

solution 

is 

Eq. 

(A14): 

(2)  When  m3 
Is  =  h  =  Bi 

>  m2  and  nx  <  n2, 

+  B2  +  ^  B3  i. 

the 

solution 

is 

Eq. 

(A14): 

(3)  When  m3 
Is  =  h=  Bi 

<  m2  and  nx  <  n2, 
+  B2-fB3i. 

the 

solution 

is 

Eq. 

(A15): 

(4)  When  Hi] 
I,  =  h  =  B, 

>  m2  and  nx  >  n2 
+  B2-^-B3i. 

the 

solution 

is 

Eq. 

(A15): 

(5)  When  m, 

h  =  Io  =  - 

=  m2  or  nx  =  n2, 

do 
*1  ' 

the 

solution 

is 

Eq. 

(A13): 

For  the  given  parameters  a,  I<,  R,  Tm,  TLX,  TH2,  TL2,  mx ,  m2,  nXl  n2, 
one  can  obtain  the  system  stable  current  7S,  and  then  obtain  the 
cooling  load  and  COP  of  the  combined  thermoelectric  device,  see 
Eqs.  (9)  and  (10). 
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